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We formulate the BES/FRS equation as a functional equation in the rapidity space and 
perform its strong coupling expansion in the limit when i = j/Ag is kept finite. We 
obtain a result which is consistent with the previous calculations at tree level and one 
loop and which differs from the two- loop calculation in 0805.4615 by a term singular at 
i = 0. 



1 Introduction and main results 



An important test of the AdS/CFT correspondence is based on the comparison of the 
anomalous dimension of the Wilson twist L operator [1] and the energy of the folded 
string spinning on AdSaxS^. The general twist L operator has the following form: 

TiiD'I'ZDI'Z ...Dl^Z), ni+n2+...+nL=M, (1) 

where is the light cone covariant derivative and Z = $i + i<p2 is a complex scalar 
field. 

In the spin chain picture the fields Z correspond to the nodes of the sl{2) spin chain 
and the covariant derivatives D to the excitations (magnons). The anomalous dimension 
of the operator ([T]) plays the role of the energy of the corresponding spin chain state. The 
energy of the ground state 7[(7,L, M] is a function only of the gauge coupling constant 
9 — ^1 the length L and the number of magnons M. This function admits a remarkable 
logarithmic scaling when M oo [1]: 

7 = /[^/,...]logM + ..., (2) 

which is valid only for certain scaling behavior of L: at weak coupling L < logM and at 
strong coupling for L < g log M. 

The important case with L/logM — » was initially considered [1,2]. In this case 
the prefactor f[g] in ([2]) depends only on the coupling constant and is equal to twice the 
cusp anomalous dimension rcuspf^f]. The cusp anomalous dimension is an object which 
can be defined for any four-dimensional gauge theory. For A/" = 4 SYM it was computed 
up to four loops [3]. 

In the limit M ^ oo, L/ logM the Bethe Ansatz equations which describe the 
sl{2) spin chain can be reduced to the linear integral equation, known as the BES equa- 
tion [2,4]. This equation allows to find the cusp anomalous dimension at any value of 
the coupling constant. The condition that the weak coupling expansion of the solution 
matches with the 4-loop calculations in the gauge theory side was an important ingredi- 
ent for establishing the complete form of the Bethe equations. Using the BES equation, 
the strong coupling expansion of the cusp anomalous dimension was performed numer- 
ically [5] and analytically at the leading [6-9] order. In [8, 10] a recursive procedure for 
analytical expansion to any desired order was given. The obtained results reproduced 
the string theory calculations at tree level [11,12], one [13] and two [14] loops. The first 
three orders of the strong coupling expansion are given by 

1 31og[2] 2^ 1 

e TT vr^ 4g 

The cusp anomalous dimension is an example of a function which smoothly interpolates 
between weak and strong coupling regimes. Freyhult, Rej and Staudacher [15] proposed 
a more general function which might had this property. They considered the limit 

M,L -.oo, J = r^ fixed (4) 
log M 
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and showed that in this case the logarithmic scahng ([2]) exists at all orders of the pertur- 
bation theory . The limit (jll) was initially introduced in [1] and the logarithmic scaling 
at g = and arbitrary j was observed in [16]. In the limit the function f^^'^[g,j] 
depends on two parameters and is called the generalized scaling function. It can be 
extracted from the solution of the linear integral equation (the BES/FRS equation), 
introduced in [15]. 

On the string theory side the logarithmic scaling is in particular realized in the 
following limiting procedure, which is taken in two steps: 

1) gr^Lr^M CXD, 

2) M/9»L/,»1. £=j^^^ fimte. (5) 

The prefactor of logM, which we will denote as f^^""-9^ jg given by the following strong 
coupling expansion 

f string ^| ^ }_^j^string ^ ^j^string ^ ^2 f string + (g) 

Comparing the limits (jlj) and ([5]), Freyhult, Rej and Staudacher raised the question 
whether the strong coupling expansion of f^^^[g,j] in the limit g,j^oo and i = j /Ag 
fixed coincides with /'^*""^[e, £]. 

From the string theory perturbative calculations the tree [11,12] and the one loop [13] 
results for the expansion ([6]) were obtained for any value of ^. At two loops only the first 
two orders of the small I expansion were found [17]: 

/f = + 4(4""' - 6 log£ + 8 log^ + 0{t) , (7) 

where got"' = - f log2 + |. 

The limit ([5]) was also used for the calculations from the asymptotic Bethe Ansatz. In 
this limit the expansion ([6]) for arbitrary value of £ was found at tree and one loop [18,19] 
and then at two loop order [20]. While tree and one loop calculations coincide with the 
string theory predictions, the small expansion of the two-loop result [20] is different 
and is given bj0: 

f2\^) = --2+ -2{<iE2^ - 6 log£ + 8 log^ I) + 0{t') (8) 

withgo^/ = -|log2 + ^. 

A different limit, particularly interesting from the string theory perspective, was 
proposed by Alday and Maldacena [16]. In this limit g ^ oo and j is exponentially 
small with respect to g: 

J ^m^ ^^/^e""^. (9) 



'^We use normalization which is different from the one used in [20]. 
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In this limit only massless excitations on are important. Therefore, the theory should 
be described by the 0(6) sigma model. The parameter m is identified with the mass gap 
of the 0(6) sigma model. In the limit the difference between the generalized scaling 
function and twice the cusp anomalous dimension can be expanded in the powers of j 
and is given by the following expression 

/^""'[g^j] - 2T,^M = -j+fE[j/m] + O(j^) . . . , (10) 

where the term j'^E[j/m] is identified with the energy density of the ground state of the 
0(6) sigma model. The corrections of the order O(j^) cannot be obtained from the 0(6) 
sigma-model. 

Basso and Korchemsky [21] applied the Alday-Maldacena limit to the BES/FRS 
equation and derived the thermodynamic Bethe Ansatz of the 0(6) sigma model at 
zero temperature. Therefore, the strong coupling expansion of the generalized scaling 
function in this limit should reproduce the string theory predictions. 

The expansion of E[j/m] at small j/m was computed in [21-23]. The expansion at 
large of j/m was done in [24] resulting in the following expression for the generalized 
scaling function: 

= -J + e'f + ^ - 2 log^) + ^ ( qg"" - 6 log£ + 8 \ogf )j+ 0{f) (11) 

withgo^/ = -|log2 + ^. 

We see that the results ([H]) and ffTT]) coincide, although they were obtained in the 
different orders of limits. 

The main goal of the current paper is to perform the strong coupling expansion of 
the generalized scaling function f^^^[g,j] with i = j/Ag fixed. The order of limits that 
we use is different from what was used for the calculation of ([7]), ([8]) and from what was 
used for the calculation of fllip . At tree level and one loop we obtain the result which 
coincides with /q*^*"^[^] and f^^"'^^[i] (and therefore with the Bethe Ansatz calculations 
in the limit ([5])). At the two-loop order we obtain the answer for arbitrary i which can 
be written in terms of f2^[(!]- 

+ (12, 

It is interesting to consider the large i expansion of /2[^]- The reason for this is the 
following. As we can conclude from [13], at the first few orders the large j expansion 
of the scaling function f[g,j] should have the BMN-like properties. This means the 
following. The large j expansion has the form 

ns.j] = E'^E"iM (13) 

n>l m>0 

and the coefficients cio, cn, cu, C20, C21 do not depend on the coupling constant. All these 
coefficients except the C12 can be found from tree [11,12] and one loop [13] calculations 
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on the string side. The prediction for them was confirmed by numerical computation at 
weak couphng [25]. In the current paper we derive the coefficient C12 from f2^^[^]- We 
find that 

fi'''[^] = lYs+0{r') and ci2 = y. (14) 

We verify this prediction by the numerical computations at weak coupling (see Appendix 
iBl) . Note that both f2^^[i] and f2^[i] give the same prediction for C12. 

The article is organized as follows. In Sec. 2 we derive the BES/FRS equation from 
the Baxter-like equation. In Sec. 3 we analyze the analytical properties of the resolvents. 
In Sec. 4 we formulate the perturbative solution at strong coupling and perform explicit 
calculations for tree level, one and two loops. For the two-loop order we use the results of 
Sec. 5, in which we analyze the behavior of the solution near the branch points. Finally, 
in Sec. 6 we summarize the obtained results. 



2 BES/FRS equation 

We start from the Bethe Ansatz equations for the sl{2) sector [4,26-28]: 




(15) 



Here e = ^ and the normalization of rapidities is suited for the strong coupling expan- 
sion. 

The variable x is the inverse Jukowsky map of u: 



1 / n . . I. /. 1 . 



u = - ^a; + -J , x[u] = M ^1 + y 1 - — I , X [u] = x[u ± ie] . (16) 

The branch of the square root is chosen in a way that |a;| > 1. 

We assume M to be even and enumerate the Bethe roots in a way that Uk > ui for 
k>l. 

The dressing phase 0[M,f] can be represented in the form 

0[u,v] = ^ {x[x^,y'] + x[x',y^] + x[x', -y'] + x[x^, -y^]) ■ (17) 

We used the notation = x[v ± ie]. We will also use y = x[v] below. 

The function x[3^5Z/] is analytic for |x| > 1 and \y\ > 1 and respects the following 
parity properties: 

x[x,y] = -x[y,x] = -x[x,~y]. (18) 

Since we consider the ground state which is symmetric, the Bethe equations will not 
change if we replace the dressing phase with the function x[a;"'",|/~] -|- x[x~ yy~^]- 
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In the following we will use the Baxter-like equation 

Q\u\ T[u] = W[u + ie]Q[u + 2ie\ + W[u - ie\Q[u - 2ie] (19) 
introduced in [19]. Here 

M 
2 



i— — in 
2 



~ 2 

VTbize] = ix^f ff fl-^^1 gT2ix[x±,-J]. (20) 

The equation (fT9|) should be understood as the definition of T[u]. The requirement of 
analyticity of T[u\ on the real axis is equivalent to the requirement for Uk to satisfy the 
Bethe equations ( ITSl) with the dressing phase 6[u,v] replaced by x[^^ ^U^] + x[^~ ^U^]- 
The equation ( JT9l) implies that T[u\ has L zeros which are usually called holes. One can 
show that both the Bethe roots Uk and the holes are real. 

The equation f|T9|) resembles the Baxter equation for the eigenvalues of the transfer 
matrix. However, there is no known transfer matrix or equivalent object with eigenvalues 
given by zeros of T[u\. 

To proceed, we introduce the resolvents 

The distribution of the Bethe roots (magnons) is given by the density function Pm[u\- 
It is supported on the two intervals [—a^xt, —<A U [a, a^xti with a ~ 1 and a^xt ~ "S*. The 
density pm[u] is related to the resolvent Rm'- 

Pm[u] = {Rm[u + ZO] -Rm[u- zO] ) . (22) 

As we argue in the appendix |Al the density is finite at u = ±a and the resolvent Rm 
has a logarithmic type singularity at these points. In other words, the typical distance d 
between the roots near u = ±a approaches zero as in the considered limit. Therefore, 
d is much smaller than the shift e in the equation (1191) for any finite value of e. This 
implies the fact that + ie\Q[u + 2ie] > W[u - ie\Q[u - 2ie] for ^[u] > and 
W[u + ie]Q[u + 2ie] < W[u - ie]Q[u - 2ie] for ^[u] < 0. 

It is instructive to compare the limit (jlj) with another well-studied limit L M oo. 
In the latter case e -C and at the distances from the real axis of order e, which are 
finite, the two terms of the r.h.s of f[T^ have the same magnitude. 

In the limit (jl]) all the holes except two are supported on the interval [—a, a]. Their 
distribution is given by the function given by the discontinuity of Rh on [—a, a\. 

Ph[u] = {Rh[u + lO] -Rh[u- 20] ) . (23) 
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There are two holes which are situated outside the interval [—aextiO'ext\- The position 
of these two external gives us the information about the asymptotic behavior of Rm for 
large absolute values of u. One can show that at the scales u ~ log M the resolvent Rm 
is constant in the leading order of the limit (jlj) [2]. In our normalization this constant 
equals to — ^. At scales larger than u ~ logM the roots do not contribute to the leading 
logM order of the energy. Therefore, we will consider the scale u ~ logM as infinity. 
The asymptotic behavior of the resolvents is given by 



R-n 



Rh 



i f3 

T- + -,u 
e u 

3 



oo ± iO; 



u 



oo. 



u 



The generalized scaling function can be found by: 

/ = -2/3-j. 



(24) 



(25) 



In the following we consider '^[u] > 0. We can neglect the second term in the r.h.s 
of ( fT9l) . Taking the logarithmic derivative of (fT9l) . we obtain the equation 



{l-D^)Rm + Rh 
where D is the shift operator 



d \ogW[u + ie] 
du log M 



0, 



(26) 



D 



(27) 



Using the resolvent Rm we can rewrite the sum over the Bethe roots as the contour 
integral. In particular, the term in fl26|) containing logiy[u + ie] can be rewritten as 



d log + ie] 
du log M 



-2D 



■^^du ( los 



R-iO 



1 

xy 



+ ix[x,y]] DRm - D 



L 



log M X du 



We used the fact that W[u + ie] as the function of Uk is analytic in the lower half plane. 

Further simplification can be achieved by performing the contour deformation, ex- 
plained in Sec. 3.2 of [8]: 



oo— iO 

dv 



2m 

'OD—iO 



<9„log 



1 - 



xy 



l+iO 



DRmlv] 



dv y ~ y 



1 



27ri X — - V 

-l+iO ^ 



u 



-DRm[v] 



■DRr, 



(29) 



where K± - the kernels introduced in [8] 



{K^F)[u] 



dv y 



^ 1 



-l+iO 



27ri X — - V 

X 



U 



{F[v + zO] ± F[~v + iOD 



(30) 
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Figure 1: Left: Analytical structure of i?m M • Right: Analytical structure of -R/i [n] . Solid line 
- the only cut on the physical sheet which corresponds to the roots/holes condensation. Dashed 
lines correspond to the cuts which appear if we analytically continue the resolvents from the 
upper half plane by the path shown by arrow. 



We can make the same contour deforming trick with duX- The conjectured BES/BHL 
dressing phase [4, 29] is such thal^l 



oo— iO 



f dv 

/ iduX[x,y]DR^[v]— = -K.^—j^K+DRm- 



(31) 



-oo— iO 



Using ([2H]), (EHD and dSU), we get from ([2SD the functional version of the BES/FRS 
equation vahd in the upper half plane: 



;i - D^)Rm + Rh = -D {K_+K+ + 2K_D 



1 \ 1 £ dx 

DkA DR^ + -D-^ (32) 



1 — D"^ J e X du 

with i = ej. This equation can be also derived from the original BES/FRS equation, 
formulated in the Fourier space [15], by means of an inverse half- Fourier transform [30]. 
The neglecting the nonlinear term in the derivation of the BES/FRS equation from the 
nonlinear integral equation [31] is equivalent to the neglecting one of the terms in the 
r.h.s of the equation f[T^ . 



3 Analytic properties of the resolvents and possible 
regimes at strong coupling 

The analytical structure of the resolvents can be deduced from (l32l) and is shown in 
Fig. [3l There are two types of the branch points of the resolvents. The first type of 
branch points has "kinematic" origin - these are branch points of the inverse Jukowsky 
map x[u] which is singular at m = ±1. As we can deduce from the equation (152]) (and 
its conjugate which is valid in the lower half plane), the "kinematic" branch points are 

^The expression j-^jyi should be understood as a series over positive powers in D. 
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located at u = ±1 + (2Z + l)ie. The branch points of the second type are located at 
u = ±a + 2Zie, where a[g,j] is the end of the root distribution. The only branch points 
on the physical sheet are u = ±a. All the others appear after the analytical continuation 
through the cut to the nonphysical sheets. 

In the following we make an assumption that the generalized scaling function is 
bianalytic function of g and j everywhere except the values of g and j for which the 
branch points of different type collide. For real g and j the collision is impossible except 
g = oo. However, if we take into account complex values of these parameters, the 
collision is possible and it determines the radius of convergence for the Taylor expansion 
of f[g,j] around some analytical point. For example, for j = we have a = 0. The 
kinematic branch point touches the origin when g = z(l/4 + Z/2). Therefore, the radius 
of convergence of the weak coupling expansion of 2rcusp[5'] equals 1/4. This observation 
coincides with the numerical prediction in [4]. 

At the strong coupling all the kinematic branch points condense onto u = ibl. There- 
fore, if we perform the strong coupling expansion of the generalized scaling function 

f[9,j] = l{fo + efi + e'f2 + ...) (33) 

the coefficients /„ should become singular when the end of the root distribution ap- 
proaches ±1. 

In the scaling limit with i = ej fixed the position of the branch point at e = is 
given by [18,20]: 

a[e = 0,e] = ^(b+^ , b = VT+i^. (34) 

Therefore, for any positive values of i the functions are analytic. And also, the 
large i expansion commutes with the strong coupling expansion. At first few orders the 
large i expansion has BMN-like properties as we mentioned in the introduction. 

At £ = the functions fn[(-] should become singular. To investigate in detail this 
special case we introduce an additional resolvent H via the relation 

(D - D-^){Rm-H) = 2K+DRm (35) 

together with the demand for Rm and H to have the same asymptotics at infinity. 
The BES/FRS equation then can be rewritten in the following forn|§ 

{D - D-'){Rm + H) = 2K^DH + 2D-^Rh - . (36) 

Let us consider the region < 1 and perform an analysis similar to one in [8]. We 
act by Ka^ on (l35l) and by on (1361) and, since K\ = K±, obtain 

= K+{{D + D-^)Rra + iD-D-')H), 

= K_ {{D-D-^)R^-{D + D-^)H -2D-^Rh). (37) 

^The equations p2|) . ((35)) and ((36|) were obtained in collaboration with Ivan Kostov and Didina 
Serban. 
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Treating this expression perturbatively, we conclude from the definition of K± and after 
some algebra that for < 1 

{D^ + D~^){Rra[u + iO] -Rm[u- ^0]) = {D- D-'){D-'Rh[u + tO] + DRh[u - ^O]). (38) 

The notion "perturbatively" means that we should perform the strong coupling expansion 
before evaluating the expression. For instance, the expression D~^Rh[u + iO] should be 
understood as 

D-^Rh[u + 20] = Rh[u + iO] - ieduRh[u + + ... (39) 

but not as D~^Rh[u + iO] = Rh[u — it + iO] = Rh[u — it]. 

In the region v? < the discontinuity of Rm is zero by definition, therefore 

D-^Rh[u + iO]+ DRh[u-iO] = . (40) 

At least perturbatively, the zero modes of D — D^^ do not contribute. One can show (see, 
for example [8]) that the zero modes of D — D^^, even if present, lead to nonperturbative 
corrections which are not considered in this papeiEl- 

For > the discontinuity of the resolvent Rh is zero. Therefore, if < < 1 
we conclude from (l38l) that 

- D-^ 

Rm[u + zO] - R^[u - zO] = ^2 _^ j^_2 ^h. (41) 

We see that the density of magnons inside the Jukowsky cut [—1, 1] has the same mag- 
nitude as the density of holes which is of order j. For a < 1 the magnitude of j is of 
order g°' e^'^3^"'~^^ [21,24] and therefore is exponentially small with respect to the coupling 
constant. We see that the perturbative expansion of the resolvent in the powers of e takes 
place only for m > 1. 



4 Perturbative solution at strong coupling 

In the following we will consider the strong coupling limit with i = ej fixed. We assume 
the following expansion of the resolvents 

Rm = -Rm,0 + Rm,l + ^Rm,2 + • • • , 

Rh = ^Rh,Q + Rh,i + eRh,2 + . . . . (42) 

We will treat all the equations in this section perturbatively in the sense that first we 
perform the strong coupling expansion of the resolvents and of the shift operator and 
then evaluate the expressions. 

■^However, they are important for the comparison with the 0(6) sigma model [21]. 
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4.1 Disappearing of the dressing phase 

As was discussed in the precedent section, in the case when C. remains finite all the roots 
are situated outside the Jukowsky cut. We will show that the BES/FRS equation can 
be considerably simplified in this case. 

We start the discussion with the following observation. Before performing the limit 
(jlj) the resolvent of the magnons can be represented as 

log[M] R,n = Y. = Z + ■ (43) 

In the continuous limit, the first term on the r.h.s will give a function which has a cut 
on the interval (—00, —b*] U [b*, 00) of the Jukowsky plane, where b* is related to a by 

The second term on the r.h.s of ( l43l) leads to the function with a cut on [—1/6*, 1/6*]. It 
is convenient to introduce the resolvent S[x] with the only branch cut (—00, —b*] U [b*, 00) 
by the relation 

Rm = S[x] + S[l/x] . (45) 

It is important to distinguish the exact position of the branch point b* with the 
position of the branch point for e = which we denote by b. 

The equality (H5l) does not fix S[x] uniquely, but up to an equivalence 

S[x] ~ S[x] + Q[x] , Q[x] + Q[l/x] = . (46) 

In particular, S[x] ^ ^ i;^^- "^^^ structure of the BES/FRS equation gives us a 

preferred choice of the representative of {S}. We demand (x^ — 1){D — D^^)S to be 
analytic everywhere except on the cut {—00, —b] U [b, 00). We will refer to this demand 
as the condition of analyticitjQ . Note that the condition of the analyticity should be 
understood perturbatively in the sense that we explained above. In particular, in the 
first three orders the condition of analyticity implies that 

5*0 is analytic outside the cut, 
5*1 is analytic outside the cut, 

5*2 d^So is analytic outside the cut. (47) 

6 

^Before performing the limit (jl]) the resolvent S[x] can be written as 

with Q[x] being meromorphic odd function with poles at x = ±1. (5[a;] is chosen to satisfy the condition 
of analyticity. 
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The preferred choice of the condition of analyticity is explained by the following pertur- 
bative equality 



K+DR^ = {D - D-')S\l/x\ (4^ 



Indeed, 

K^DR^ = f^ ^'l ^ iD-D''){S[y] + S[l/y]) 
J Zm X — - V — u 



1 1 



-1 



^^^(D - D--)S[l/y] = {D- D-^)S\l/x]. (49) 

X 

The contour of integration goes clockwise around [—1,1]. We used the condition of 
analyticity to close the contour of integration. 
From (HHj) and (l35l) we deduce that 

H = S[x]-S[l/x]. (50) 

The BES/FRS equation simplifies: 

(1 - D^)S + Rh = -DK_{D - D-^)S + -D-^. (51) 

e X du 

This simplification is quite remarkable. The form of flSTl) is similar to the BES/FRS 
equation ( 132|) with DRm substituted with {D — D^^)S and significant simplification of 
the kernel which becomes just K^. We see that the convolution of the kernels, which 
comes from the term with the dressing phase in the Bethe Ansatz, has disappeared at 
strong coupling. This is an all-loop generalization of the observation, made in [20], that 
the strong coupling expansion of the Bethe equations does not contain S'[l/x] terms. 

The equation (ISTl) itself cannot be solved since it contains two unknown quantities. 
Together with it we should use the conjugate equation, which is valid in the lower half 
plane: 

1 / rlr 

(1 - D-^)S + Rh = -D-^K(D ~ D-')S + -D'^-^. (52) 

e X du 

If we multiply ( ISTj) by -D^^, fl52|) by D and take the difference of the obtained equations 
on the interval > 6^, we will obtain the following result^ 

S[x + iQ] + S[x-iQ] = -Rh. (53) 

Let us discuss the relation of the equations (!5T]) and (!53l) with the BES equation, which 
we expect to recover when Rh = and i = 0. The equation (15T1) is a variation of the idea 
to decouple the BES equation into two simpler equations by introducing an additional 

^Strictly speaking, we obtain {D — D^^){S[x + iO] + S[x — iO] + Rh) = 0. Asymptotics at infinity 
assures that it is equivalent to (fSSjl . 
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function [8,10]. From the equalities fH5|) and fl50|) we can identify S with the resolvents 
in [8]: 

S[l/x] = R+, S[x] = R^. (54) 

This identification is possible only perturbatively at strong coupling. We see that at 
strong coupling the resolvents -R+ and R- are not independent but can be expressed 
through the one resolvent S, the physical meaning of which we discussed above. 
For Rh = the equation fl53|) simplifies to 

S[x + iO] + S[x - iO] = 0. (55) 

Its solution in the leading order reproduces the AABEK density [7]. This equation in 
another parametrization was discussed in [8]. 

Of course, the equation (1551) has an infinite number of solutions. The correct solution 
is chosen from the investigation of the analytical structure of 5* in the NFS double scaling 
limit [8,10]. 

Note that we obtained the equation S[x + iO] + S[x — iO] = in a different order of 
limits than what was used in the strong coupling expansion of the BES equation: we 
took i after the strong coupling expansion while the BES equation corresponds to 
the inverse case. Although the order of limits do not affect the equation itself, the 
analysis in the vicinity of the branch point is completely different in two cases. When 
we apply the NFS double scaling limit we automatically use the order of limits for the 
strong coupling expansion of the BES equation. 

4.2 The perturbative procedure 

Using the equations fl5T]) . fl52|) . and fl53l) . we can perform the strong coupling expansion 
of the resolvent S. The logic of the computations is the following: since D — D^^ = 0{e), 
using ( 15T|) and ( l52l) . we can express Rh,n in terms of Sm with m < n. Then, by solving 
the Riemann-Hilbert problem (!53|) for Sn-, we can express Sn in terms of Rh,n- The final 
result is given by the equations (162!) which allows to express Sn in terms of Sm with 
m < n. 

In the following we discuss the details of the calculations which lead to (1^2]) . 
The equation (IHTi) allows us to find Rh[u] in the upper half plane. Let us introduce 
the function $ by the following relation 

$ = D-^Rh - %u] > 0. (56) 

e X du 

From flSTl) we conclude that we can rewrite the function $ in the following way 

l+iO ^ 

^x] = iD-D-')S[x]- I ^ ^'\ \ {D-D-'){S[y]-S[l/y]) = 

J iTXl X — - V — U 
-1+iO 

- (O - , / (-1- - -ij) (Z, - O-^Sm . (57) 
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where F = r_|_ + r_, r± - the upper/lower unit semicircle from —1 to 1. 

We perform the deformation of the contour of integration, shrinking it towards the 
interval [—1,1]. Since '^[u] > 0, the contour r_ passes through the point l/x['u]. The 
residue at this point cancels the first term in the r.h.s of (!57|) and we obtain 



dy y - y\ I 1 



2m \x — - I \y — X y 



{D-D-')S[l/y] 



7 = 7+ + 7-, 7± = [-l±«0,l±zO] . 



(5^ 



We take (1581) as the definition of $ for any values of x. From the equation ( ]52l) we get 
that 



1 / dr 
e X du 



(59) 



As it follows from ( l58ll . the function is discontinuous everywhere on the real axis of 
the Jukowsky plane. However, from ( l59l) and ( l56i) we can expect that 



D-^<^[x -iQ] = D^x + x^>h'^ 



(60) 



since Rh should be continuous on the interval x"^ > b"^. Indeed, we can prove the equality 
(l60l) by induction using the definition fl58l) of $. 

Using (I56l) . fl59|) . and fl60l) we can reformulate fl53|) as an equation for the function S: 

S[x + iO] + S[x - iO] = ^[x + iO] + $[a: - iO] + , x^ > b\ (61) 



The general solution of ( ]6Tll which respects the cut structure of 5* is the following: 



Sn — Snh,n + Sho 



1 2 

7l> + l>-i 



2n+l 



om.n 
Snh,n 



h,n + h,n 
1 



_ 1 ' 

X 

Cn,k 



(1 _^2)fc + (52 _ ^2)fc 

fc=l ^ ^ fc=l ^ ^ / 



D + D-^ I 2Tri\x-^ 



i] {D-D-')S[l/y]x 



1 



y - X 
1 



^52 _ 3,2 



/ 1\ 2 

dy y ~ y\ 1 



n-1/2 



D + D-^J 27ri\x-y y - i (62 - 



.(i^-Z?-i)5[l/y] 



• (62) 



Here /3 = /3+ + /?_ with /3± = [-1/6 ± iO, 1/b ± iO]. 
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We remind that we assume the perturbative expansion of S 



S = - {So + eSi + e^S2 + ...)■ 



(63) 



By the notation [F]n we understand the coefficient in front of e"~^ of the perturbative 
expansion of F in powers of e. The perturbative expansion of the integrals Ji and I2 
is understood in the following way: first we expand the resolvent and the expressions 
containing the shift operator D and then perform the integration. Since D — D~^ = 0{e), 
the r.h.s of (|62l) contains only Sm with n < m. Therefore, the solution (1621) defines the 
recursive procedure which allows to express Sn in terms of Sm with n < m. 

The constants an,k are fixed by the condition of analyticity, which in the first three 
orders is given by fHTl) . 

To fix the coefficients Cn^k we need an additional information. At the ffist two orders 
these coefficients are fixed by the following properties: 

• Sn has at most a — x"^)^^"^ " singularity at x = ±6 (in particular, 5*1 [±6] = 0), 

• 5*1 behaves as — ^ for x 00 + iO, 

• Sn decreases at infinity for n > 1. 

The asymptotics at infinity fixes c„^i = — a„.i for any n. 

These conditions are not sufficient to fix the coefficients Cn,k with n > 2 and A; > 1 
which appear starting with two loops. The analysis which allows to fix these coefficients 
is given in Sec. [51 

The generalized scaling function, determined by the formula fl25|) . can be expressed 
in terms of the coefficients an,k and Cn^k by 



The asymptotics at infinity fixes ao,i = — 1 and the condition of analyticity is satisfied 
only if 




(64) 



In the following we will focus on the ffist three orders. 



4.3 Tree and one loop level 

The leading order solution follows easily from 0621) : 




(65) 



X 



(66) 



We obtain 



(67) 
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At the one loop level the solution fl62|) reads 



Si 

Shom,l 
Snh,l 



Snh,l + SjiomAi 



X 



x^ 



+ 



Oil 



IlA + /2,1, 



X — - 



dy_ 

2-K : 



(1 - X2)2 62 _ 3,2 



-^5o[l/l/]J X 



1 



i2 — ^2 



^2,1 



y - X 

2 

T 

X 



VP 



27r 2 



-|;S„|1M 




\ _ a;2 



(6^ 



The integrals /i^i and /2,i can be evaluated explicitly (see appendix [C] for the details) 
and the coefficients ai^i, ai,2 are fixed by the condition of analyticity of 5*1 at x = ±1. 
The resolvent 5*1 is given by the following expression: 

Snh,i = -2l{\og[b - 1] - 21og[&] + log[6 + 1]) + 

e 



+ 



(41og[6] - 21og[6 + 1] - \og[b^ + 1]) + 



+ l[x+ - ]\og 



+ 



+ 



{I - x)(^/W^^ + I x) 

(1 +x)(\/62 -x2 X) 

(1 -x){h + -x2 + ^ x) 



+ 



\/62 — x2 



62 



log 



log 



{I + x){h - x){VW — x^ — £ x) 



+ 



52 



62 



(69) 



ai,i 



ai,2 



(36^ _ 352 + 5 _ 1 _ 452 (^^2 _ 3^ ^Qg[^] ^ (52 _ 3) log [6 + 1] + 



+ [h" - 1) log [6 - 1] + - (6^ _ 452 _ 1) iog[^2 ^ 



1 



(4-46-8 log[6] + 4 Iog[6 + 1] + 2 log[62 + 1]) 



The one-loop generalized scaling function is given by 

6-1 + 862 log[6] - (62 + 1) iog[52 + 1] - 2(6^ - 1) log[6 - 1] - 26^ log[6 + 1] 



A 



27r 6 



(70) 



(71) 



The results for /o[£] and fi[P\ coincides with the results in [18,19] and therefore with 
the perturbative calculations in the string theory. The difference between 5*1 and the 
resolvent found by Belitsky [19] is an odd meromorphic function of x with poles at 
X = ±1. This difference does not contribute to Rm- 
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4.4 Two loops 

The two loop resolvent is given by the solution ( 162|) with n = 2. The computation at 
this order is more difficult since we have to fix the coefficient C2,2 in Shom,2 which requires 
an additional analysis of the initial equation (132!) . We postpone this analysis to the next 
section. The result is the following: the resolvent 5*2 is singular at the point x = b with 
the leading square root singularity given by 

Q 

= (26)3/2 (5 _ a;)3/2 + • • • C^^) 



with 



e = ^(-2 + 26 - 66^ + 6b^ - 46^ log[2] - 46^ (26^ + 6 + 2) log[6] + 

+2 (6^ - 26^ - 26^ + 1) log[6 - 1] + (46^ + 46^ + 46^) log[6 + 1] + 

+ (6^ + 462-l)log[62 + l]). (73) 

This singularity comes from Shom,2 and Snh,2- The contribution from these two terms is 
the following 

^hom,2 (26)3/2 (6 - a;)3/2 
6^ P 

Snh,2 = (26)3/2 (6_a;)3/2 +•••' C^^) 

where P is a number given by a complicated integral. We give the explicit form of P in 
the appendix [Dl 

Comparing (!72|) with (17i|) we get the following expression for the constant 02,2: 

C2,2 = b\Q-P). (75) 

The form of the solution suggests that in the vicinity of x = 1 the function Snh,2 can be 
expanded as 

Snk,2 = + + 0{1). (76) 

(x — 1)"' (x — 1) 

The coefficients 02,^ are fixed from the condition of analyticity which implies that 

+ (^^ + (^ " 2^ + 3^7 ^^^^ 

should be analytic at x = 1. 
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Collecting all the coefficients together, we get the following expression for the 2-loop 
generalized scaling function 

Exphcit expression for —2ilsi + 4£/s2 (l + 2I2) is given in the appendix [Dl 

The numerical comparison with the result of [20] shows that the scaling function fl75]) 
and the scaling function obtained from the Bethe Ansatz in the order of limits ([5]) are 
related as 

fr'[i] = f2^[(]+m, (79) 

11 11' 



m 



ViTP V4:i^ 12 £4 



The expression ( 1781) can be expanded for the large values of i. We are interested in 
the leading term from which we can find the coefficient C12 defined in (JT3l) . Evaluating 
first 



3 TT^ TT 

32 _ Jlog[2] + log[£])log[2£]^ ^ ]_ ^ 

-2iu + UL,{i + ^)f = (-6 + ^)^ + (^(r^), (80) 

we obtain 

/2M = i^ + 0(r^) and ci2 = y. (81) 

The term 5[^] does not contribute to the coefficient C12. To verify ( IHTl) . we performed 
numerically large j expansion of the generalized scaling function at (7 = 0. The details 
of the computation can be found in the appendix [Bl 



5 Behavior of the solution near the branch point. 

In this section we explain how to obtain the coefficient Q in fl72|) . 

The perturbative solution ( !62|) is defined for the values of x such that \x — h\ is much 
larger than e. At the point h the perturbative expansion is not valid and the solution 
( l62l) develops a singularity. Let us understand what type of the singularity is expected 
at this point. For this we take the difference of (15T!) and (!52|) for x'^ > b"^, gettingj 

1 2/ 

D S[x + iO] + D-^S[x - iO] =K_{D- D-^)S ^. (82) 



^Strictly speaking, we get {D — D ^) (|82|) = 0. The equation ((82|) is deduced from the conditions of 
decreasing at infinity 
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Figure 2: The scales that appear in the Bethe equations, b is the position of the smallest 
positive Bethe root, d is the distance between two neighboring roots, e appears as the 
shift in (HM. 



The r.h.s of this equation is analytic at x = b and therefore we can write the following 
equality: 

D Ssing[x + tO] + D-^Ssing[x - tO] = 0, (83) 

where by Ssing we denoted a part of the resolvent 5* which is singular ai x = b. We 
will be only interested in the leading singularities of S at each order of the perturbation 
theory. At tree and one loop level they can be found directly from the corresponding 
solutions fl^S]) and flB^ . We find the most general form of the two loop singularity by 
solving the equation (1831) . Combining the leading tree, one and two loop singularities 
together, we get the following expression: 

_ Vb- X \fW' { log[6 - x] + 



1 (tb^\ ^ Q - 4(log[6 - x] + 0) + (log[6 -x] + 0)2 



2\7rey {b-x] 



2 



+ . . . 



= (-2 + 26 - 66^ + 66^ - 46^ log[2] - 46^ (26^ + 6 + 2) log[6] + 

+2 (6^ - 2b^ - 2b^ + 1) log[6 - 1] + (46^ + Ab^ + 46^) log[6 + 1] + 

+ (54 + 462-i)log[62 + i]). (84) 

The coefficient Q is arbitrary. To fix Q we need to consider the BES/FRS equation in 
the double scaling regime which is defined as follows. We introduce the variable 

. = ^ (85) 
e 

and perform the double scaling limit e — ^ with z fixed. Note that a is the exact position 
of the branch point. 

Since e ^ (i, the treatment of the resolvent as the analytic function of z with cuts 
is valid in the double scaling limit. To compare, if we apply the double scaling limit for 
the order of limits ([5]), used for the perturbative calculations in the string theory, we will 
see the separate poles of the resolvent since e <^ d. We illustrate the discussion in this 
paragraph in Fig. 2. 
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For our purposes it is sufficient to keep only the leading order of the double scaling 
limit. From the solution ( l62i) we see that the leading order scales as e"^''^. The functions 
analytic at the point u = a scale with integer powers of e and therefore will not appear 
in the double scaling limit. 

We see that the equations (l5Ti) and (l52l) reduce to the following ones 

S\z\ - S\z + 2i] + Rh[z\ = foru.h.p, 

S[z] - S[z - 2t] + Rh[z] = forl.h.p. (86) 

To solve these equations we perform the inverse Laplace transform. We define the inverse 
Laplace transform for S and Rh in the following way: 

— ioo— 



S[s]= [ ^S[z]e^\s<0, 
J 2m 



ioo— 
— too+O 

~ f dz 

Rh[s]= I —Rh[z]e'%s>0. (87) 

— ioo+O 

We chose different contours of integration for S and Rh in order to avoid the cuts of 
the resolvents. The functions S[s] and Rh[s] are defined for any complex value of s by 
analytical continuation of fISTI) . The cut of S[z] implies the existence of the cut for S[s] 
which we define to be on the ray s > 0. Equivalently, there is a cut for Rh[s] which we 
define to be on the ray s < 0. 

The inverse Laplace transform of the equations fl86|) gives the following equations 
valid for s > 0: 

(1 - e-^'')S[s - tO] + Rh[s] = 0, 

{1 - e+^'')S[s + iO] + Rh[s] = 0. (88) 
They imply the equation on S[s] 

S[s + iO] = -e-'^''S[s - tO], s>0, (89) 

which is solved by 

Here A[s] is a function with trivial monodromy. The factor r[l — -]e~^ can be absorbed 
in A[s]. The reason to keep it explicitly will be clear below. 

At large z the function S'[z,e] should be glued with the solution (1^ . In particular, 
the leading e ^ order of S[z, e] in the double scaling limit should reproduce fl8^ . This 
demand is satisfied if A[s] is expanded in positive integer powers of s around s = 0: 

A[s] = a{l + ais + a2s'^ + . . .). (91) 
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Indeed, assuming flOTl) we get from the Laplace transform of S[s] the following large z 
expansion for S[z]: 

s[z] = f f: % +f: +f: ^^-^^"5'^^' + • • -V 02) 

^ \ Z Z Z I 

\ n=0 n=l n=2 / 

The coefficients do,™ are expressed in terms of with m < n and dfc^n are expressed in 
terms of (io,m with m < n — k. Upon identificatioi 

a = ^=^, 6* = 6 + -(e-21og2-log7r) + 0(e') (93) 



2&2 e ' v^£3' ^ 

the expression (!92l) coincides with (18^ with 

Q = TT^ - a2 + 2a2^ . (94) 

We used the substitution z = -^^-^ instead of flHSl) which is possible for the consider- 
ation of only the leading singularities. The overall normalization a is chosen to fit ([22]) 
with the singularity of 5*0. The exact position of the branch point b* is chosen to fit (!92|) 
with the singularity of ^i. The coefficients ai and 02 are still arbitrary. They define the 
unknown coefficient Q through flM|) . 

To fix «! and 02 we use the fact that S[z] has a logarithmic singularity at the origin. 
Therefore, S[z] is expandecj^ at the origin as 

S[z] = (ro + nz + r2Z^ + . . .) log[-2;] . (95) 

This type of the singularity comes from the logarithmic singularity of the resolvent at 
finite coupling. Another way to see this singularity is to take the difference of the first 
and the second equations in ( l86ll for 2; > 0. The resulting equation is 



S[z + iO]-S[z- iO] = S[z + 2i]~S[z- 2i] . (96) 

We see that the discontinuity of S[z] is given by the function S'[2; + 2i] — S[z — 2i] which 
is analytic at z = 0. This implies the expansion (p5!) for the function S[z]. 

The series (19 5 p has finite radius of convergence due to the branch points z = 2iL of 
S\z\ in the nonphysical plane. This means that S\s\ admits only an asymptotic expansion 
at infinity 

5^W = -7 + 5-25 + --- (97) 

which is valid in the cone d < arg s < 27r — -(9 with arbitrarily small 'd. In the same cone 
the asymptotic expansion of 74[s]/S'[s] can be performed 

.e^--'°g[-^l = — 1 + — - + —— + ... . (98) 



r[l - f ] ^2 V 12 s 288 s 



^We remind that b* is the exact position of the branch point of the resolvent S[x] and b = VT+P. 
^the expansion of the type ro + n.^ + r2z'^ may contribute as well, but only in subleading orders. 
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Therefore, A[s] in this cone also has an asymptotic expansion 



A[s]=qo + ^ + ^ + .... (99) 



s s"^ 



On the other hand, using (IHHl) . we can express A[s\ in terms of Rh'- 

A[s] = -— r [-1 s^/^e^~^'°^^^^Rh[s]. (100) 

Repeating the same arguments we conclude that A[s] admits asymptotic expansion (p9|) 
in the cone — (tt — < args <'K — d with arbitrarily small d. Therefore, the series (!99|) 
is valid everywhere providing with this to be convergent. 

Since R[s\ is analytic outside s < and S[s\ is analytic outside s > 0, from the 
comparison of fllOOp and ( |90i) w e deduce that A[s\ is analytic in C*. 

This allows us to identifjPl ([99]) with We conclude that A[s\ is a constant and 
therefore 

Q=\t,\ (101) 
From (18^ and fllOip we get the expression (1731) for Q. 



6 Conclusions 

In this paper we computed the strong coupling expansion of the generalized scaling 
function form the BES/FRS equation up to two loops. At the two-loop level the result 
is different from what was obtained in [20]. A possible reason for the discrepancy is in 
the different order of limits that was used in two approaches. The BES/FRS equation 
is derived in the limit (jlj) in which we can neglect one of the terms in the l.h.s of the 
equation ( JT9l) . On the other hand, in the order of limits ([5]) used in [20], this is not 
justified near the branch point of the resolvent. The order of limits is certainly not 
important for the calculations at the tree and the one-loop level since the analysis in the 
vicinity of the branch point is needed starting from the two loops. 

We also gave a prediction for the leading term of the large £ expansion of f2[C] flHT]) 
and checked it numerically using the BMN-like properties of the expansion f[T^ . In 
view of the discrepancy between f2^[i] and f2^""'^[i] it would be interesting if one can 
reproduce the result (18T1) from string theory calculations. 
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A The structure of the resolvent in the vicinity of 
the branch point 

Our aim is to show that the resolvent in the vicinity of the branch point has the loga- 
rithmic cut. For the sake of simplicity we discuss the Bethe eqautions for g = 0. We 
choose the normalization of u such that the Bethe equations at zero coupling constant 
have the form 

Uk - y fji Uk -Uj+t 

If we take the logarithm of the Bethe equations, we obtain 

M 

Uk - Uj] = sign [uk] , (103) 

i=i 

. , 1 u + i 

F[u\ = - — : log ■ 



27ri u — i 

The equation fll03p can be interpreted as the force equihbrium equation in classical 
mechanics. For the distances between the particles much larger than one, the interaction 
between the particles can be approximated with the Coulomb force F[uk — Uj] ~ ^ uk-u- • 
In this case the density of the particles is approximated by the square root cut in the 
leading order of the large spin limit. 

However, for the scaling considered in the current paper we have the opposite situa- 
tion: the distances between the particles with the smallest absolute values of rapidities 
are much smaller than one. Therefore, the Coulomb approximation is not applicable. To 
describe the distribution of roots in the considered limit it is better to represent F[u] as 

F[u] = arctan[M] H — sign [u] . (104) 

TT 2 

M 

If we introduce the effective force -Feff[-u] = L F[2u] — J2 't arctan['u — Uj], then for the 

i=i 

positive roots equation of the equilibrium will be written as 

i^eff['"Af/2+fe] = ^ + ^ (105) 

Since is a smooth function, we immediately get that in the vicinity of the branch 
point 

The last estimation comes from the dominant LF[2m] term of the and is consistent 
with assumption of the small distance between the roots. It is valid for Uuj^ ^ which 
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is the case for any finite j. Equidistant distribution between the roots corresponds to 
the logarithmic type of the branch point in the continuous limit. 

The situation does not change for any finite value of the coupling constant. The 
additional terms in the Bethe Ansatz do not change the arguments used in derivation 
because the interaction introduced by them is nonlocal for any finite g. 

At the infinite value of the coupling constant two effects appear. First, the roots 
scale generically as g and therefore the distance between them becomes large. In fact, 
in the order of limits ([5]) used for the perturbative calculations in the string theory, the 
branch point changes into a square root type for any nonzero i. 

Second, at strong coupling and for i the distribution of roots approaches the 
Jukowsky branch point, where the additional terms of the Bethe equation become local. 
They result in the change of square root cut into the u~^^^ behavior of the resolvent. 



B Large j expansion 

Large j expansion for the Bethe equations f|T5l) at (7 = was done numerically in the 
regime (jlj) [25]. However, only the first two terms of this expansion were given. We need 
the third term in order to verify the prediction for the coefficient C12 = 

The generalized scaling function can be found if we know the density of holes ( 
[15,25]): 



f{g,j) = Sg'' + 2g' I duU[l + iu]+ij[^-iu] - 2^(1) ) + 



(107) 



The density of holes satisfies the integral equatio: 



Ph=l-j:^(i^[l + iu]+ij[^-iu]j+j^{ij[i{u-v)]+^[-i{u-v)^^ (108) 

—a 

and is normalized by 

J= [ Ph{v)dv. (109) 



Substituting the normalization condition into f llOSp we get an integral equation which 
depends only on the parameter a. Solving numerically this equation and fitting the 
results for j in the range from 30 to 150, we obtain the following large j expansion: 



f{g,3) = g\— — + ^-...) + o{g^), (110) 



,8.0000 6.79 5.33 

j 3 

which is consistent with ( ISTj) . 

^^to the 5 = Bethe equation 

^^it is derived under the same assumptions as the BES/FRS equation. 
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C Evaluation of the Cauchy type integrals 



The solution (1621) suggests evaluation of the Cauchy type integrals. This evaluation can 
be simplified if we will use the observation explained below. 

Suppose that we know the value of the following Cauchy type integral 



I\x] 



dz 



[III] 



z — X 



with 7 - any contour, closed or not. Then for any rational function Q[x\, regular on the 
contour 



F[x\ = dz 

J ■y 



z — X 



Q[x]I[x] +R[x] 



(112) 



with R[x] being the rational function. 

Indeed, discontinuous part of F is given by. 

F[x + iO] - F[x - iO] = Q[x]f[x] = Q[x]{I[x + iO] - I[x - iO]). 

The most general solution of the last equation is 

F[x] = Q[x]I[x] + R[x] 



(113) 



(114) 



with R[x] being a rational function. Since F[x] should be the value of the integral (11121) . 
it should be analytic outside the contour 7 and decrease at infinity. R[x] is chosen to 
fulfill this properties. 

We use it to compute integrals in (1^ . The building blocks (i.e. integrals of the type 
I[x]) for calculation of 1 loop correction can be just guessed after some experience and 
are given by 



dy- 



1 



y-x 



log 



X — 1 

x+l 



1/6 



dy- 



-i/h 
1 



/ 



dy- 



y-x 
1 



log 



X — 

X + 



dy- 



y-x 
1 



VP 



x^ 



dy 



y-x 
1 



x^ 



x^ 



\J\p- — x^ log[ 



X 



x^ 



-xV^^+Vb^-x^' 



^ log[ 



x^ 



62 



1 

-71 



log I 



x^ 


jb^ - 




-Vb^ 


- 1 
— 1 


X y 


jb^- 


1 

X'' 


+ Vb' 


- 1 


x^ 


jb^ - 


1 


- Vb^ 


— 


X \ 


h - 


1 


+ Vb^ 


— x"^ 



(115) 
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where L = [-1,-1/6] U [1/6,1]. 

One can check that the discontinuity of the r.h.s on the contour of integration coin- 
cides with the integrand in the l.h.s. 

In view of the relation (11121) . knowledge of the integrals (11151) makes the calculation 
of Si straightforward. 



D Integrals for the 2-loop correction 

Direct calculations shows that 



7+ 



/3+ 



p = P1 + P2, 

£ f dy jy' - l)ib' - j.r^'hiSl[l/y]) - Bo 
' 27r 1-6V 

0+ 

Bo = 4^(464a;log[- ] + {-2 + 2h - Qb' - 2^^ + 2{l - 2h'' - 2}P + h'^)\og[b - I] - 

2l^TT OX +1 

862(1 + h") log[6] + (4fe2 + Ah^ + 4fe4) log[l + 6] + (-1 + 46^ + 6^) i^gfi + 52^^^^ 
^2 = ^^^^(-2 + 26-662-263 _63log[4] + 2(1- 2&2-2fe3 + 64)log[6-l] - 

26^(4 + 6 + 46^) log[6] + (46^ + 46=^ + 46^) log[l + 6] + (-1 + 46^ + 6"*) log[l + 6^]). (117) 
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